Heat transport by Dirac fermions in normal/superconducting graphene junctions 
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We study heat transport in normal/superconducting graphene junctions. We find that while the thermal con- 
ductance displays the usual exponential dependence on temperature, reflecting the s-wave symmetry of the 
superconductor, it exhibits an unusual oscillatory dependence on the potential height or the length of the bar- 
rier region. This oscillatory dependence stems from the emergent low-energy relativistic nature of fermions in 
graphene, essentially different from the result in conventional normal metal/superconductor junctions. 



The recent progress in practical fabrication techniques for a 
monoatomic layer of graphite, called graphene, has allowed 
for experimental studies of this system, which in turn has 
triggered a tremendous interestii^iii^i^. Graphene is a two- 
dimensional system of carbon atoms, and the low-energy elec- 
trons in graphene are governed by Dirac equation. Up to 
now, intensive studies on graphene have been conducted for 
instance quantum Hall effect*'*^''^, minimum conductivity^ and 
bipolar supercurrentiS. 

From applied physics point of view, graphene is also an 
important material. Graphene exhibits high mobility and car- 
rier density controllable by gate voltage, which makes it well 
suited for achieving device applications . ^•^•"•'^ In order to ap- 
ply graphene to electric devices, it is an important issue to 
clarify characteristics of transport phenomena in graphene. 

In conventional normal metal/superconductor junctions, it 
is known that electric and thermal conductances reflect the 
magnitude or symmetry of the gap of the superconductorjJ^iii 
While conductance in normal/superconductor graphene junc- 
tion has been studied, '^''^''^ thermal conductance in the same 
junction has not yet been investigated. The study of the ther- 
mal conductance in normal/superconductor graphene junction 
will complement the study of the conductance in the same 
junction. 

In this paper, we study heat transport in nor- 
mal/superconducting graphene junctions. We find that 
the thermal conductance has an exponential dependence on 
temperature, which reflects the s-wave symmetry of the su- 
perconductor. However, it displays an oscillatory dependence 
on the potential height or the length of the barrier region. 
This oscillatory dependence stems from the relativistic nature 
of fermions in graphene, and differs in an essential way from 
the result in the conventional normal metal/superconductor 
junctions. 

We briefly present the formalism to be used in this pa- 
per, following Ref.'^. Consider a two dimensional nor- 
mal/insulating/superconducting graphene junction'** where 
the superconducting (normal) region is located in the semi- 
infinite regions a; > L (x < 0). The proposed experimental 
setup of our model is shown in Fig. [T] By exploiting the val- 
ley degeneracy,'^ the Bogoliubov-de Gennes equation for the 
junction in the xy-plane reads 
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FIG. 1: (color online) The proposed experimental setup to 
measure heat transport by Dirac fermions in a graphene nor- 
mal/superconductor proximity structure. The top and bottom gate 
allow for the chemical potential in the middle region to be adjusted. 
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The superconducting order pa- 
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where Q{x) is the Heaviside step function, while ip is the 
phase corresponding to the globally broken U{1) symmetry 
in the superconductor Also, vp ~ lO^m/s is the energy- 
independent Fermi velocity for graphene, ai{i — x,y) de- 
notes the Pauli matrices, E is the excitation energy, and u and 
V denote the electron-like and hole-like exictations, respec- 
tively, described by the wave-function. The Pauli matrices op- 
erate on the two triangular sublattice space of the honeycomb 
structure, corresponding to the A and B atoms. The linear dis- 
persion relation is a reasonable approximation even for Fermi 
levels as high as 1 eV,^" such that the fermions in graphene be- 
have like massless Dirac fermions in the low-energy regime. 

Let us consider an incident electron from the normal side 
of the junction {x < 0) with energy E. For positive excitation 
energies _E > 0, the eigenvectors and corresponding momen- 
tum of the particles read 



r+ = [l,e'«,0,0]^e 



p'' = {E + Ef)/vf, (3) 



for a right-moving electron at angle of incidence 9, while 
a left-moving electron is described by the substitution 6 
TT — 6. The superscript ^ denotes the transpose. If Andreev- 
reflection takes place, a left-moving hole with energy E and 
angle of reflection 6a is generated with corresponding wave- 
function 



(1) 



V'- = [0,0,1, 
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{E-Ef)/vf,(4) 
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where the superscript e (h) denotes an electron-Hke (hole- 
like) excitation. Since translational invariance in the y- 
direction holds, the corresponding component of momentum 
is conserved. This condition allows for determination of the 
Andreev-reflection angle 9a via p'' sin 0a ~ p'^sinO. From 
this equation, one infers that there is no Andreev-reflection 
{9a = ±7r/2) for angles of incidence above the critical angle 

9c^8in-WE-EF\/{E + EF)). (5) 

On the superconducting side of the system {x > L), the pos- 
sible wavefunctions for transmission of a right-moving quasi- 
particle with a given excitation energy E > reads 

^gig=cose + ^^ ^ (^/^ ^ y/E^-A^)/vF, (6) 
^gig'-cose-x^ ^ (^/^ _ y/E^-A^)/vF. (7) 

The coherence factors are given by^i 

Above, we have defined 61+ = 9g and6'~ = """^s- Thetrans- 

mission angles 6*^ for the electron-like and hole-like quasi- 

particles are given by g^'-' sin 9^"^ = p'^ sin 6, i=e,h. Note that 
in all the wavefunctions listed above, for clarity we have not 
included a common phase factor e''^"^ which corresponds to 
the conserved momentum in the y-direction. 



It is appropriate to insert the restriction which will be used 
throughout the paper, namely A <C Ep. Since we are using 
a mean-field approach to describe the superconducting part 
of the Hamiltonian, phase-fluctuations of the order parameter 
have to be small^\ 

We define the wavefunctions in the normal, insulating and 
superconducting regions by ip, ipj and ^t, respectively, with 

V; = V4+rV'l+?-AV'-, (10) 

i^i = hr+ + 12^- + + (11) 

* = (12) 

The wavefunctions ip differ from V' in that the Fermi energy 
is shifted by an external potential, such that Ef — > Ep — U 
where U is the barrier height. Also, note that the trajectories 
of the quasiparticles in the insulating region, defined by the 
angles 9 and 9a, differ by the same substitution: 

sin9/sm9 ={E + Ef)/{E + Ef-U), (13) 
sm9A/sin9={E + EF)/{E-EF + U). (14) 

Note that the subscript ± on the wavefunctions indicates the 
direction of momentum, which is in general different from the 
group velocity direction. 

By matching the wavefunctions at both interfaces, V'U=o = 
'0/U=o and "0/12:=^ — '^\x=L we obtain the following 
expressions for the normal reflection coefficient r and the 
Andreev-reflection coefficient r^jii 

r^teiA + C)+thiB + D)-l, (15) 
rA = te{A' + C') + thiB' + D'), (16) 

where the transmission coefficients read 



te = 2cos6l[e-"^-*(B' + D') - {E'c''^^ - D'c'^^)]p-\ (17) 
th = te[e'^^{A'e''^^ - C"e'^^) - A' - C"][B' + D' - e'"-^ (B'e"'^"-^ - D'e'^'^)r\ (18) 
p = [e-'^^(B' + D') - {B'e-'^^ - D'e'^-')][e-'\A + C) + {Ae'^ - Ce"'^^)] 
- [{De-'^ - - c-'%B + D)] [A'c"'^^-^ - C'c'^-' - c'"'^^ [A' + C')] (19) 

I 

and we have introduced the auxiliary quantities 

A = we'(«^-P^'[l - (e'^^ -e'^^)(2cos^)-i], 
B = we'^'^-P^^il - (e'^~ - e"'")(2cos0)-i], 
C = we'(P^+'?^)(e'*^ - e'^^)(2cos0)-\ 
D = i;e'(P^+«")(e'^^ - e'^")(2cos^)-\ 



A' = z;e'(«*+f"-*)[l + {e'^* -e-'^^){2cos9AY\ 
B' = we'(9"+P"~'^)[l + (e'*" - e-'^^){2coB9Ay\ 
C" = z;e'(«*-P"-*)(e-'^^^ - e'^^)(2cos^A)"\ 
D' =ue'^^'-^''-'l'\e-'^'' -e'^")(2cos^A)"^ (21) 



(20) 
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Here, we have defined 

(7+ = COS0+L, = q'^ cosO^ L, 
= cos 6L^ ~ p^ cos Oa^- 



with X = LU/vp. This indicates that thermal conductance is 
TT-periodic with respect to x in this limit. 



(22) 



In the thin-barrier limit defined as L ^ and U ^ oo, one 
gets 



0, Oa -^0, q±^ 0,p± TX 



(23) Finally, the normaUzed thermal conductance is given by 



J-tt/2 



(24) 



r 



with the gap at zero temperature Aq = A(0). 





1 ' 1 ' 1 

^kFL=0 


' / 




----kFL=10 






- ■ ■ -kFL=100 






1 -r^' i h- 


(a) ; 


H 


^kFL=0 


— 1 — ^ 




-----kFL=10 






■ - kFL=100 








(b) . 

1 



0.2 0.4 0.6 O.S 



0.6 
0.5 
0.4 
5^ 0.3 
0.2 
0.1 





n ' r- 



r, ■/■-■-- 
\ ' . i-' 
\ i 
■ (a) 



-kpL=0 
-kpL=l 
kpL=2 
-kpL=5 



2 




4 6 


8 10 




1 ,_ 






- - - -U/Ep=0 
• ■ ■ • •U/Ep=l 

U/Ep=2 

U/Ep=5 



10 



20 



FIG. 2: (Color online) Thermal conductance as a function of T/Tc 
for various kpL with U/Ef = 10 and E'p = lOOAo at Ep = 
lOOAo in (a) and Ep = lOAo in (b). 



FIG. 3: (Color online) (a) Thermal conductance as a function of 
U/Ep for various kpL with T/Tc = 0.5 and Ep = E'p ^ lOOAo. 
(b) Thermal conductance as a function of kpL for various U/Ep 
with T/Tc = 0.5 and Ep = E'p = lOOAo. 



We next present our results for the normalized thermal 
conductance. Figure |2] (a) shows thermal conductance as a 
function of T/Tc for various kpL with U/Ep = 10 and 
Ep = E'p = IOOAq. Here Tc is the transition temper- 
ature and hp = Ep/vp. From Fig. |2] (a), an exponen- 
tial dependence of the thermal conductance on temperature is 
seen, similar to the conventional normal metal/superconductor 
junctions.'^ This exponential dependence reflects the s-wave 



symmetry of the superconductor However, the length depen- 
dence of the thermal conductance is nonmonotonic (oscilla- 
tory) and thus essentially different from that in the conven- 
tional normal metal/superconductor junctions. A similar plot 
for Ep = IOAq is shown in Fig. |2](b). We also find an ex- 
ponential temperature dependence, but the dependence on L 
gets weaker Therefore, the magnitude of the oscillation with 
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respect to kpL gets reduced with the increase of the Fermi 
wave vector mismatch. 

Figure |3] (a) depicts thermal conductance as a function 
of U/Ep for vai'ious length kpL with T/Tc = 0.5 and 
Ep = E'p = lOOAo. An oscillatory dependence of the 
thermal conductance on U/ Ep is seen. The period decreases 
with kpL. Figure [3] (b) displays thermal conductances as a 
function of kpL for various U/Ep with T/Tq = 0.5 and 
Ep = E'p — IOOAq. We also find an oscillatory dependence 
on kpL. The period also decreases with These fea- 

tures stem from the 7r-periodicity of the thermal conductance 
with respect to x = kpLU/Ep in the limit of U ^ Ep 
and kpL <^ 1 similar to the junction conductance i '^'^^ In 
other words, the damped oscillatory behavior of the thermal 
conductance is a direct manifestation of the relativistic low- 
energy Dirac fermions. Also, the presence of the insulating 
region is essential for the oscillatory behavior. 

Since we have assumed a homogeneous chemical poten- 
tial in each of graphene regions, the experimental observa- 
tion of the predicted effects require charge homogeneity of 
the graphene samples. This is a challenge, since electron- 
hole puddles in graphene imaged by a scanning single electron 
transistor device^^ suggest that such charge inhomogeneities 
play an important role in limiting the transport characteris- 
tics of graphene^*. In addition, we have neglected the spatial 



variation of the superconducting gap near the interface. The 
suppression of the order parameter near the interface is ex- 
pected to be least pronounced when the sharp edge criteria is 
satisfied and there is a large Fermi-vector mismatch. In the 
present case, this is precisely so, whence we do not expect 
our qualitative results to be affected by taking into account 
the reduction of the gap near the interface. Finally, we have 
assumed that there is no lattice mismatch at the interfaces and 
that these are smooth and impurity-free'*. A more refined pic- 
ture could be obtained by using more realistic models of the 
variation of the chemical potential, i.e. a continuous slope in- 
stead of a step-like variation. 

In summary, we have studied heat transport in nor- 
mal/superconducting graphene junctions. We found that the 
thermal conductance has an exponential dependence on tem- 
perature which reflects the s-wave symmetry of the supercon- 
ductor but oscillatory dependence on the potential height or 
the length of the barrier region. This oscillatory dependence 
stems from the relativistic nature of fermions in graphene, es- 
sentially different from the result in the conventional normal 
metal/superconductor junctions . 
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